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Abstract 



o 
l> 

o 

^ I We study the average number of simplices N'{r) at geodesic dis- 

tance r in the dynamical triangulation model of euclidean quantum 
gravity in four dimensions. We use N'{r) to explore definitions of cur- 
Oh! vature and of effective global dimension. An effective curvature Ry 

goes from negative values for low K2 (the inverse bare Newton con- 
stant) to slightly positive values around the transition k^. Far above 
the transition Ry is hard to compute. This Ry depends on the dis- 

K> ' tance scale involved and we therefore investigate a similar explicitly 

C^ ■ r dependent 'running' curvature Res{r). This increases from values 

of order Ry at intermediate distances to very high values at short 
distances. A global dimension d goes from high values in the region 
with low K2 to d = 2 at high H2- At the transition d is consistent with 
4. We present evidence for scaling of N'(r) and introduce a scaling 
dimension ds which turns out to be approximately 4 in both weak 
and strong coupling regions. We discuss possible implications of the 
results, the emergence of classical euclidean spacetime and a possible 
'triviality' of the theory. 
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1 Introduction 

The dynamical triangulation model ll|, |],|[ is a very interesting candidate 
for a nonperturbative formulation of four- dimensional euclidean quantum 
gravity. The configurations in the model are obtained by glueing together 
equilateral four- dimensional simplices in all possible ways such that a sim- 
plicial manifold is obtained. A formulation using hypercubes was pioneered 
in 0]. The simplicial model with spherical topology is defined as a sum over 
triangulations T with the topology of the hypersphere S"^ where all the edges 
have the same length i. The partition function of this model is 

Z{N,^,)= Yl exp(«:2iV2). (1) 

T,N4=N 

Here iVj is the number of simplices of dimension i in the triangulation T. The 
weight exp(fi;2A''2) is part of the Regge form of the Einstein- Hilbert action 
-J^R/IQttGo, 

-S = j^Y.V226a = K2{N,-pN,), (2) 

.2 = ^, ,^10^^e^o<V4)^ 

8Go 271 ' ^ ^ 

Here V2 is the volume of 2-simplices (triangles A) and ^a is the deficit angle 
around A. The volume of an z-simplex is 



v, = r^it + i)/2yt\. (4) 

Because the Ni {i = ... 4) satisfy three constraints only two of them are 
independent. We have chosen N2 and A^4 as the independent variables. For 
comparison with other work we remark that if A'^o is chosen instead of A'"2 
then the corresponding coupling constant Kq is related to ^2 by kq = 2fi;2. 
This follows from the relations between the Ni, which imply that 

No - ^N2 + N, = x, (5) 

where x is the Euler number of the manifold, which is 2 for S"^. 



Average values corresponding to ([1|) can be estimated by Monte Carlo 
methods, which require varying A^4 [|l],^|^. One way to implement the con- 
dition iV4 = A^ is to base the simulation on the partition function |]I|] 

Z'{N, t,2) = E ^M^2N2 - ^4N, - 'y{N, - N)'), (6) 

r 

where k,4 is chosen such that (^4) ^ N and the parameter 7 controls the 
volume fluctuations. The precise values of these parameters are irrelevant if 
the desired A4 are picked from the ensemble described by (|^). This is not 
done in practice but the results are insensitive to reasonable variations in 
7. We have chosen the parameter 7 to be 5 ■ 10~^, giving (A^|) — (^4)^ ~ 
(27)"^ = 1000, i.e. the fluctuations in A4 are approximately 30. 

Numerical simulations [^0,0,^,|^,|^,|^ have shown that the system (^ can 
be in two phasescl. For K2 > K2{N4) (weak bare coupling Go) the system is in 
an elongated phase with high (R), where (R) is the average Regge curvature 

In this phase the system has relatively large baby universes [Q and resembles 
a branched polymer. For K2 < '^2(^4) (strong coupling) the system is in a 
crumpled phase with low (R). This phase is highly connected, i.e. the average 
number of simplices around a point is very large. The transition between the 
phases appears to be continuous. 

As an example and for later reference we show in figure |1| the susceptibility 
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(R) - (R)' V. (8) 



{V = A4V4). The two curves are for A^ = 8000 and 16000 simplices. 

The data for A^ = 8000 are consistent with those published in reference 
where results are given for 



x{No) = ^ ((A^o') - {No)') . (9) 



*In section we discuss the possibility that these are not phases in the sense of con- 
ventional statistical mechanics. 
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Figure 1: The curvature susceptibility as a function of K2 for 8000 and 16000 
simplices. 
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For fixed A^4 this is 1/4 of our curvature susceptibihty (|^), as can be seen 
from (I). 

The behavior of Z{k,2,N) as a function of A^ for large A^ has been the 
subject of recent investigations |1^, |TT|, |13|, |12[ . In ref. [0 we discussed the 
possibility that ^2 might move to infinity as A^ ^ 00 and argued that this 
need not invalidate the model. So far a finite limit is favoured by the data 
[0,|l2|, however. 

It is of course desirable to get a good understanding of the properties of 
the euchdean spacetimes described by the probability distribution exp(— S"). 
A very interesting aspect is the proliferation of baby universes 0. Here we 
study more classical aspects like curvature and dimension, extending previous 
work in this direction ||1|,^|^,||,||,|^,^. The basic observable for this purpose 
is the average number of simplices at a given geodesic distance from the 
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Table 1: Number of configurations used in the various calculations. 

arbitrary origin, N'{r). We want to see if this quantity can be characterized, 
approximately, by classical properties like curvature and dimension, and if 
for suitable bare couplings K2 there is a regime of distances where the volume- 
distance relation N'{r) can be given a classical interpretation. It is of course 
crucial for such a continuum interpretation of N\r) that it scales in an 
appropriate way. 

In section ^ we investigate the properties of an effective curvature for a 
distance scale that is large compared to the basic unit but small compared to 
global distances. Effective dimensions for global distances are the subject of 
section ^ and scaling is investigated in section ^ We summarize our results 
in section ^ and discuss the possible implications in section ^. 



2 Curvature 

A straightforward measure of the average local curvature is the bare curvature 



at the triangles, i.e. (R) = A8Tr^j3/5i~'^{N2/Ni-p) (this follows from (0) and 
(^). It has been well established that at the phase transition {N2/N4 — p) ^ 
2.38 - 2.10 ^0.28, practically independent of the volume [|l],|,|,|]. This 
means that (R) ~ 55£~^ has to be divergent in the continuum limit £ — *■ 0. 
However, the curvature at scales large compared to the lattice distance i is 



not necessarily related to the average curvature at the triangles. One could 
imagine e.g. a spacetime with highly curved baby universes which is flat at 
large scales. 

The scalar curvature at a point is related to the volume of a small hyper- 
sphere around that point. Expanding the volume in terms of the radius of 
the hypersphere results in the relation 

V{r) = a.r-{l--^^ + 0{r% (10) 

^" = r(n/2 + l)' ^^^^ 

for an n-dimensional manifold. We have written Ry here to distinguish it 
from the small scale curvature at the triangles. Differentiating with respect 
to r results in the volume of a shell at distance r with width dr, 

V'{r) = C„nr"-i(l - ^^ + 0{r^)). (12) 

on 

We explore this definition of curvature as follows. We take the dimension 
n = 4, assuming that there is no need for a fractional dimension differing from 
4 at small scales. For r we take the geodesic distance between the simplices, 
that is the lowest number of hops from four-simplex to neighbour needed to 
get from one four-simplex to the other. Setting the distance between the 
centers of neighbouring simplices to 1 corresponds to taking a fixed edge 
length in the simplicial complex of vTO, i.e. we will use lattice units with 
i = v/lO. For V{r) and V'{r) we take 

V{r) = V,sN{r), V'{r) = V^sN'{r), N'{r) = N{r) - N{r - 1), (13) 

where N[r) is the average number of four-simplices within distance r from 
the (arbitrary) origin, N'{r) is the number of simplices at distance r and we 
have allowed for an effective volume Ves per simplex which is different from 
V4. Since Ry is to be a long distance (in lattice units) observable we shall 
call it the effective curvature. 

The effective curvature was determined by fitting the function N'{r) to 

N'{r) = ar^ + br^. (14) 
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Figure 2: The number of simplices N'{r) at distance r from the origin at 
K2 = 0.80, 1.22 and 1.50, for A^ = 16000. 



It then follows from ([T2|) and (|T3]) that Ry and Ves are determined by 

Rv = -2A-, Kff=^. (15) 



The constant C4 in equation ( PUJ ) is 7r^/2. In fiat space we would have 
Ves = V^, giving 

a = -^ = ^^ - 8.47, (16) 

where we used V4 = 25-\/5/24 for £ = V^- Such a space cannot be formed 
from equilateral simplices because we cannot fit an integer number of sim- 
plices in an angle of 2tt around a triangle. We expect therefore Ves to be 
different from V4. 

Figure 1^ shows N'{r) for 16000 simplices. Three different values of k,2 are 
shown, 0.8 (in the crumpled phase), 1.22 (close to the transition) and 1.5 (in 
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Figure 3: Effective curvature fit and global dimension fit in the crumpled 
phase at K2 = 0.80. 



the elongated phase). These curve can also be interpreted as the probability 
distribution of the geodesic length between two simplices. Such distributions 
were previously presented in refs. [§|,|§]. 

We determined N'{r) by first averaging this value per configuration suc- 
cessively using each simplex as the origin. We then used a jackknife method, 
leaving out one configuration each time, to determine the error in Ry- 

Our results in this paper are based on A^ = 8000 and 16000 simplices. 
Configurations were recorded every 10000 sweeps, where a sweep is defined 
as N accepted moves. The time before the first configuration was recorded 
was also 10000 sweeps. We estimated the autocorrelation time in the average 
distance between two simplices to be roughly 2000 sweeps for A^ = 16000 and 
K,2 = 1.22. For other values of k,2 the autocorrelation time was lower. The 
number of configurations at the values of K2 used are shown in table |^. 




Figure 4: Effective curvature fit and global dimension fit near the transition, 
K2 = 1.22. 



Figures ^|-g show effective curvature fits (continuous lines) in the crum- 
pled phase {n2 = 0.80), near the transition {^2 = 1.22) and in the elongated 
phase (k2 = 1.50), together with global dimension fits (see the next section) 
at longer distances, for A^ = 16000. The curves are extended beyond the fit- 
ted data range, r = 1-11, to indicate their region of validity. A least squares 
fit was used, which suppresses the lattice artefact region where N'{r) is small 
because it is sensitive to absolute errors rather than relative errors. 

For K2 < ^2 t^^ fi^s were good even beyond the range of r used to deter- 
mine the fit (except obviously when this range already included all the points 
up to the maximum of N'{r)). This can be seen in figure |^, where the fit is 
good up to r = 15. The fit with ar^ + br^ does not appear sensible anymore 
for K,2 values somewhat larger than the critical value K2, because then N'{r) 
goes roughly linear with r down to small distances. This can be seen quite 
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Figure 5: Effective curvature fit and global dimension fit in the elongated 
phase at ^2 = 1.50 (the effective curvature fit is not appropriate here). 

clearly in figure || (so the Ry fit in this figure should be ignored). 

Figure ^ shows the resulting effective curvature Ry as a function of n2 for 
8000 and 16000 simplices. For N = 8000 the fitting range was r = 1-9. We 
see that Ry starts negative and then rises with K2, going through zero. In 
contrast, the Regge curvature at the triangles {R) is positive for all K2 values 
in the figure. 

The value of a in (|15]) varied from about 0.9 at k,2 = 0.8 to 0.4 near 
the transition. These numbers are much smaller than the fiat value of 8.47 
in equation ([T6|), indicating an effective volume per simplex Ves ~ 20 — 50, 
much larger than V4 ~ 2.3. This is at least partly due to the way we measure 
distances. The distances are measured using paths which can only go along 
the dual lattice and will therefore be larger than the shortest paths through 
the simplicial complex. 
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Figure 6: The effective curvature Ry as a function of K2 for 8000 and 16000 
simplices. 

Tliere is a strong systematic dependence of Ry on tlie range of r used 
in the fit. For example, for the A^ = 16000 data in figure |^ we used a least 
squares fit in the range r = 1-11. If the range is changed to r = 1-9 the data 
for Ry have to be multiplied with a factor of about 1.5. We can enhance this 
effect by reducing the fitting range to only two r values and thereby obtain a 
'running curvature' Res{r) at distance r. We write N'{r) = a{r)r^ + b{r)r^, 
N'{r + 1) =a{r){r + lY + b{r){r + lf and define i?eff(?" + |) = -246(r)/a(r), 
which gives 



i?eff(r + i)=24 



(r + 1)3 - r^N'{r + l)/N'{r) 
{r + If ~ r^N'{r + 1)/N'{r)' 



(17) 



Figure |^ shows the behavior of -Rcff(^) for various k,2- It drops rapidly from 
large values ~ 4.5 (which is near (R) ?« 5.5) at r = to small values at 
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Figure 7: The effective curvature -Reflf(^) for various K2 and 16000 simplices. 



r ^ 8. For K2 < ^2 ^^e curves have a minimum and around this minimum 
the values of Res average approximately to the Ry^s displayed earher. 

The idea of Rcs{r) and Ry is to measure curvature from the correlation 
function N'{r) by comparing it with the classical volume-distance relation 
for distances r going to zero, as long as there is reasonable indication for 
classical behavior at these distances. Clearly, we cannot let r go to zero 
all the way because of the huge increase of Rcs- This seems to indicate a 
'planckian regime' where classical behavior breaks down. 



3 Dimension 

One of the interesting observables in the model is the dimension at large 
scales. A common way to define a fractal dimension is by studying the 
behaviour of the volume within a distance r and identifying the dimension d 
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Figure 
of r for 



3: ln[{N{r)/N{r — l))]/ln[r/(r — 1)] < > dlnV/dlnr as a function 

some values of k,2. 



if the volume behaves like 

V{r) = const, x r'^. (18) 

Such a measurement has been done in ||1],0|, using the geodesic distance and 



a distance defined in terms of the massive scalar propagator [|1J]. Arguments 
against the necessity of such use of the massive propagator were raised in 



ref. |]T6[. Although we feel that the issue is not yet settled, we shall use here 
the geodesic distance as in the previous section. 

If the volume does go like a power of r, the quantity 



d 



dlnV 
dlnr 



In N{r) -lnA^(r-l) 



(19) 



ln(r) — ln(r — 1) 

would be a constant. Figure H shows this quantity for some values of «2 
for a system with 16000 simplices. We see a sharp rise at distances r = 
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1, 2 independent of K2 which is presumably a lattice artefact. The curves 
then continue to rise until a maximum where we may read off an effective 
dimension, which is clearly different in the crumpled phase {n2 = 0.80-1.20) 
and in the elongated phase (^2 = 1.30-1.50). Instead of a local maximum 
one would of course hke a plateau of values where d In V/d In r is constant and 
may be identified with the dimension d. Only for K2 beyond the transition a 
range of r exists where dlnV / dhir looks like a plateau. In this range, d ^ 2. 

Similar studies have been carried out in 2D dynamical triangulation 
where it was found that plateaus only appear to develop for very large num- 
bers of triangles |]I5|. Our 4D systems are presumably much too small for 
dlnV/dlnr to shed light on a fractal dimension at large scales, if it exists. 
However, we feel that the approximate plateau in the elongated phase with 
d = 2 should be taken seriously. 

As we are studying a system with the topology of the sphere S"^, it seems 
reasonable to look whether it behaves like a d- dimensional sphere S'^. For 
such a hypersphere with radius tq, the volume behaves like 

V'{r) = dC,ry {sm-Y-\ (20) 

This prompts us to explore (^) as a definition of the dimension d, we shall 
call it the global dimension. For small r/ro this reduces to (|T8|). On the 
other hand for d = n the form ( pOD is compatible with the effective curvature 
form (0), with 

«v - ^. (21) 

To determine the dimension d we fit the data for N'{r) to a function of 
the form (pO|), 

N'{r) =c {sin— Y-\ (22) 

The free parameters are tq, d and the multiplicative constant c. It is a priori 
not clear which distances we need to use to make the fit. At distances well 
below the maximum of N'{r) (cf. figure ^) the effective curvature fit appears 
to give a reasonable description of the data, but it will be interesting to 
try ( pOD also for these distances. Small distances are of course affected by 
the discretization. This is most pronounced at low K2 where the range of 
r-values is relatively small. On the other hand, for small K2 the fits turn out 
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to be good up to the largest distances, indicating a close resemblance to a 
hypersphere, while at larger K2 the values of N'{r) are asymmetric around 
the peak (cf. fig. ^ and fits turn out to be good only up to values of r not 
much larger than this peak. The system behaves like a hypersphere up to the 
distance where N'{r) has its maximum, which would correspond to halfway 
the maximum distance for a real hypersphere. Above that distance it starts 
to deviate, except for small K2 where N'{r) is more symmetric around the 
peak and the likeness remains. 

For K,2 = 0.8 (figure |^) the global dimension fit ( p2|) was performed to 



the data at r = 7-21, for K2 = 1.22 (figure g) to r = 4-20 and for K2 = 1.50 
(figure 1^) to r = 3-60. 

The two descriptions, effective curvature at lower distances and effective 
dimension at intermediate and larger distances, appear compatible. Notice 
that at K,2 = 0.8 in the crumpled phase the local effective curvature Ry is 
negative while the global structure resembles closely a (positive curvature) 
sphere with radius tq = 7.6 (ro = 2rm/7r, with Vm the value where N'{r) 
is maximal). At K2 = 1.22 near the transition the effective curvature and 
effective dimension descriptions appear to coincide. At ^2 = 1.50, deep in 
the elongated phase, the effective curvature fit does not make sense any- 
more, its r-region of validity has apparently shrunk to order 1 or less. The 
effective dimension fit on the other hand is still good in this phase and the 
power behavior ([T8|) with d ^ 2 is extended by (p2|) to intermediate distances 
including the maximum of N'{r). 

Figure |^ shows the global dimension as a function of k,2 for the total 
volumes of 8000 and 16000 simplices. For small values of K2 it is high and 
increases with larger volumes. For values of K2 beyond the transition it 
quickly goes to two, confirming the statement made earlier that in this region 
N'{r) is approximately linear with r down to small r. 

A most interesting value of the dimension is the one at the phase tran- 
sition. To determine the value of Kg where the transition takes place we 
look at the curvature susceptibility of the system, figure |l|. For 8000 sim- 
plices the peak in the susceptibility is between k,2 equals 1.17 and 1.18 where 
the dimensions we measured are 4.2(1) and 3.8(1). For 16000 simplices the 
peak is between 1.22 and 1.23 where the dimensions are 4.2(1) and 3.6(1). 
Therefore the dimension at the transition is consistent with 4. As can be 
seen from these numbers the largest uncertainty in the dimension is due to 
the uncertainty in Kg. The effective dimensions have some uncertainty due 
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Figure 9: The global dimension as a function of K2 for 8000 and 16000 sim- 
plices. 

to the ambiguity of the range of r used for the fit. Near the transition this 
generates an extra error of approximately 0.1. 

4 Scaling 

To get a glimpse of continuum behavior it is essential to find scaling behavior 
in the system. We found a behavior like a d dimensional hypersphere for r 
values up to the value r^ where N'{r) is maximal. This suggests scaling in 
the form 



N'{r)=rt-'f{ — ,d) 



(23) 



i.e. N'{r) depends on k,2 and A^ through d = d{K,2,N) and r™, = rm{K,2,N). 
The occurence of d in this formula is unattractive, however, since it is ob- 
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Figure 10: The scaling function p for H2 = 0.8 at A^ = 8000 and ^2 = 1-0 at 
A^ = 16000. 
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tained by comparing N'{r) to sin'^~^(r/ro) at intermediate scales, which is 
a somewhat imprecise concept. We would like a model independent test of 
scaling. 

Consider the probability for two simplices to have a geodesic distance r, 



p{r) 



N'{r) 

N 



, 1 = ^p(r) ^ / drp{i 

r=l '' 



(24) 



which depends parametrically on ^2 and N . It seems natural to assume 
scaling for this function in the form 



1 T 

p(r) = — p{ — , r), / dxp{x, r) = 1, 



(25) 



where r is a parameter playing the role of d in ( psf ) which labels the different 
functions p obtained this way. For instance, r could be the value pm = p(l, t) 
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Figure 11: The scaling function p for K2 = 1.18 at A^ = 8000 and K2 = 1.22 
at A^ = 16000. 



at the maximum of p. This may give give problems if pm does not change 
appreciably with K2 (similar to d in the elongated phase). Other possibilities 
are r = fp{rm) or r = r'^/r'' for some k with r'^ = ^^p{r)r^. In practice 
we may also simply take t = n2 — ^2(^) ^^ some standard choice of N and 
compare the probability functions with the p(r) at this A^. 

Matched pairs of p{x, r) for A^ = 8000 and 16000 are shown in figures p!0| - 
|l^, respectively far in the crumpled phase, near the transition and in the 
elongated phase. For clarity we have left out part of the errors. Scaling 
appears to hold even for ^2 values we considered far away from the transition. 

The values of k,2{N) of the matched pairs in figs. pi]|-p^ are increasing with 
A^ in the crumpled phase and decreasing with A^ in the elongated phase. This 
suggests convergence from both sides to K2(^) ^^ ^ increases. For current 
system sizes K2(^) i^ ^^i^^ very much dependent on A^, a power extrapolation 
estimate [§] gives K2(^^) ~ 1-45. 
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Figure 12: The scaling function p for ^2 = 1-5 at iV = 8000 and ^2 = 1-3 at 
A^ = 16000. 



We can use the matched pairs of p to define a scahng dimension dg by 



^^ = «C% 



(26) 



where a and ds depend only on r. Using non-lattice units, replacing the 
integer r^ by rm,/(£/vl0) where r^ is now momentarily dimensionful, we 
can interpret (^6]) as 

N^i"-^) , (27) 



which shows that the scaling dimension characterizes the dimensionality of 
the system at small scales £ — ;> with r^ fixed. Returning to lattice units, 
taking for r^ the integer value of r where N'{r) is maximal and using an 
assumed error of 0.5, these scaling dimensions would be 4.5(3), 3.8(2) and 
4.0(1) respectively for figures |1^, [Tl| and 0. The largest errors in these figures 
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probably arises due to the uncertainty in the values of K2 we need to take to 
get matching curves, i.e. to get the same value of r. As we do not have data 
for a continuous range of n2 values, we have to make do with what seems 
to match best from the values we do have. Nevertheless, the values of dg 
far away from the transition are strikingly close to 4 when compared to the 



values of the global dimension d, which are 7.8 and 2.0 for figures |T^ and |T2. 

The scaling form ( p3l) is in general incompatible with (^), except for 

dg = d. The evidence for dg = 4: points instead to a scaling behavior of the 

form 



N'{r)=rif{-,T), (28) 

'm 

with /(x, r) = a{T)p{x, r). This further suggests scaling of the volume V'{r) 
at distance r with an effective volume Vcs per simplex (cf. ([TB|)) depending 
only on r. 

A precise definition of the scaling form p(x, r) may be given by 

p{x, t) = rmp{rmx), K2 = K2{N) such that rmpir-m) =T, N ^ 00, (29) 

where we used r = pm for illustration. Intuitively one would expect the con- 
vergence to the scaling limit (^) to be non-uniform, with the large x region 
converging first, and there may be physical aspects to such nonuniformity. 
The scaling analogue of running curvature (|T7p is given by 

~ 2 T^ / 24 3 — (i In p/(i In X , , 

Rix) ^ r^ R.,{xr^) = - ^^J^^^^JJ^^- (30) 

Fig. shows this function for the matched pair of fig. [l^ in the transition 
region. We have also included the curve for K2 = 1.23 at A^ = 16000. The 
curves do not match in the region around x = 0.5 and Res is apparently a 



sensitive quantity for scaling tests. Still, fig. |T^ suggests reasonable matching 
for a value of ^2(16000) somewhere between 1.22 and 1.23 and even the steep 
rise as far as shown appears to be scaling approximately, with i?( 16000) 
somewhat below /?(8000). We find similar scaling behavior for the matched 
pair in the elongated phase. The number of our K2 values in the crumpled 
phase is too limited to be able to draw a conclusion there. 

The steep rise appears to move to the left for increasing A^ (a scaling 
violation). A most interesting question is whether the onset of the rise (e.g. 
the X value where R = 50) continues to move towards a; = as A^ — »• 00. Such 
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Figure 13: Scaling form r^Resixr^) versus r/r^ near the the transition, for 
A^ = 8000, K2 = 1.18 (middle) and A^ = 16000, K2 = 1.22, 1.23 (lower and 
upper). The hypothetical limiting form corresponding to S^ is also shown. 



behavior is needed for a classical region to open up from x around 1 towards 
the origin x = 0. In other words, the 'planckian regime' would have to shrink 
in units of the size r^ of the 'universe'. Then r'^Ry could be defined as a 
limiting value of -R, x — *> 0. 

Since we are looking for classical behavior in the transition region it is 
instructive to compare with the classical form of R corresponding to the 
sphere 5*^, for which p = Pm sin^ 9,9 = 7rx/2, 



R 



187r2 tan9-9 



9^ 5tan9-39 



37r2(l 



ISvr^ 
T20 



x" + 



(31) 



This function is also shown in fig. O. Our current data are evidently still far 
from the hypothetical classical limiting form (|3T|). 
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5 Summary 

The average number of siniplices N'{r) at geodesic distance r gives us some 
basic information on the ensemble of euchdean spacetimes described by 
the partition function (|l]). The function N'{r) is maximal at r = r^ and 
f'mN'{r)/N shows scaling when plotted as a function of r/rm- We explored a 
classical definition of curvature in the small to intermediate distance regime 
based on spacetime dimension four, the effective curvature Ry- We also ex- 
plored a description at global distances by comparing N'{r) with a sphere of 
effective dimension d. Judged by eye, the effective curvature fits and effec- 
tive dimension fits give a reasonable description of N'{r) in an appropriate 
distance regime (figs. ||-§D. The resulting Ry depends strongly on the fitting 
range, which led us to an explicitly distance dependent quantity, the 'run- 
ning' curvature Res{r). This dropped rapidly from lattice values of order of 
the Regge curvature at r = to scaling values of order of the Ry found in 
the effective curvature fits at intermediate distances. A preliminary analysis 
of scaling behavior then suggested the possibility of a classical regime with a 
precise definition of r'^Ry in the limit of large N . We shall now summarize 
the results further, keeping in mind the ambiguity in Ry as derived from the 
fits to N'ir). 

For small K2 the effective curvature is negative. Furthermore the system 
resembles a d-sphere with very large dimension d which increases with the 
volume. This suggests that no matter how large the volumes we use, there 
will never be a region of r where the power law V{r-) oc r'^ gives good fits 
over large ranges of r. In other words, the curves for d\n.V / dXnr in figure ^ 
will never have a plateau. This behavior is consistent with that of a space 
with constant negative curvature, where the volume rises exponentially with 
the geodesic distance for distances larger than the curvature radius and if 
we look at large enough scales the intrinsic fractal dimension equals infinity. 
The resulting euclidean spacetime cannot be completely described as a space 
with constant negative curvature as such a space with topology S^ does not 
exist, and finite size effects take over at still larger distances. 

At the transition the spacetime resembles a 4-dimensional sphere with 
small positive effective curvature, up to intermediate distances. 

For large K2 the system has dimension 2. In this region it appears to 
behave like a branched polymer, which has an intrinsic fractal dimension 



of 2 0]. Moving away from the transition, the curvature changes much 
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more rapidly than in the small K2 phase and the effective curvature radius 
rv soon becomes of the order of the lattice distance. A priori, two outcomes 
seem plausible. In the first, the system collapses and ry becomes of order 1 
in lattice units, reflecting the unboundedness of the continuum action from 
below. In the second, the spacetime remains 4-dimensional and ry can still be 
tuned to large values in lattice units, but very small compared to the global 
size Tfn, for a sufficiently large system. So far the second outcome seems to be 
favored, for two reasons. Firstly, the function N'{r) looks convex for r < 6 
slightly above the transition, e.g. for k,2 = 1-24 at 16000 simplices. In other 
words, the linear behavior as seen in figure ^ does not set in immediately 
above the transition. Secondly, the system shows scaling and the scaling 
dimension dg as defined in (^B|) is approximately 4 even far into the elongated 
phase, indicating 4-dimensional behavior at small scales. 

High and low dimensions in the crumpled and elongated phases with the 
value four at the transition were reported earlier in ref. 0. This dimension 
was apparently interpreted as a small scale dimension, whereas instead we 
find a small scale dimension of four in all phases. Similar results are also 
found in the Regge calculus approach to quantum gravity ||1^|18|, where one 



also finds two phases, a strong (bare Newton) coupling phase with negative 
curvature and fractal dimension four, and (using an R^ term in the action 
for stabilization) a weak coupling phase with fractal dimension around two. 

6 Discussion 

The evidence for scaling indicates continuum behavior. This brings up a 
number of issues which need to be addressed in a physical interpetation of 
the model. 

One guideline in this work is the question whether there is a regime of dis- 
tances where V'{r) = V^fjA^'(r) behaves classically for suitable bare Newton 
constant Gq. The connection with classical spacetime can be strengthened 
by identifying the geodesic distance r with a cosmic time t and V'{ry^^ with 
the scale factor a{t) in a euclidean Robertson- Walker metric. The classical 
action for a{r) is given by 



fdr[a(^) +a] + \(27i^ fdra^-vY (32) 



2 
^~ 8G 
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where A is a Lagrange multiplier enforcing a total spacetime volume V. It 
plays the role of a cosmological constant which is just right for getting volume 
V. For positive GX the solution of the equations of motion following from 
(|32D is a = rosinr/ro, which represents S"^ with R = 12/rQ = 192'7rG'A. 
For negative GX the solution is a = tq sinh t/tq which represents a space of 
constant negative curvature R = —12/rQ = 1927rG'A, cut off at a maximal 
radius to get total volume V. 

The form (^) serves as a crude effective action for the system for in- 
termediate distances around the maximum in V'{r) at r^, and couplings 
K2 oc Gq^ in the transition region. At larger distances the fluctuations of the 
spacetimes grow, causing large baby universes and branching, and averaging 
over these may be the reason for the asymmetric shape of V'{r). Because 
of this the Robertson- Walker metric cannot give a good description at these 
distances. 

Intuitively one expects also strong deviations of classical behavior at dis- 
tances of order of the Planck length a/G, assuming that a Planck length 



exists in the model. A proposal for measuring it was put forward in [|T£ 
The steep rise in the running curvature Re«{r) at smaller r indeed suggests 
such a planckian regime. It extends to rather large r but it appears to shrink 
compared to r^ as the lattice distance decreases, i.e. N increases for a given 
scaling curve labelled by r. This suggests that the Planck length goes to 
zero with the lattice spacing, G/r^ — >• as iV — i> oo at fixed r. This does 
not necessarily mean that the Planck is length of order of the lattice spacing, 
although this is of course quite possible. However, the theory may also scale 
at planckian distances and belong to a universality class. It might then be 
'trivial'. 

At this point it is instructive to recall other notorious models with a 
dimensionful coupling as in Einstein gravity, the 4D nonlinear sigma models. 
The lattice models have been well studied, in particular the 0(4) model for 
low energy pion physics (see for example [0). It has one free parameter 
K = i'^f^ which corresponds to the renormalized dimensionful coupling /^; 
/ is the pion decay constant or the electroweak scale in the application to 
the Standard Model. With this one bare parameter it is possible to tune 
two quantities, f /m and if, where m is the mass of the sigma particle or 
the Higgs particle. This trick is possible because the precise value of i is 
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unimportant, as long as it is sufficiently smalP. However, in the continuum 
limit if ^ triviality takes its toll: m/ f -^ and the model becomes 
noninteracting. The analogy /^ <-> l/G, m ^^ r~^ suggests that we may be 
lucky and there is a scaling region in K2-N space, for a given scaling curve 
(given r), where the theory has universal properties and where we can tune 
G/r^ to a whole range of desired values. Taking the scaling limit however 
might lead to a trivial theory with G/r"^ = 0. In case this scenario fails 
it is of course possible to introduce more parameters, e.g. as in R^ gravity, 
to get more freedom in the value of G/r"^. This then raises the question of 
universality at the Planck scale. 

We really would like to replace r^ by Ry in the reasoning in the previous 
paragraph, since we view Ry as the local classical curvature, provided that 
a classical regime indeed develops as A^ — > 00. 

Next we discuss the nature of the elongated phase. Even deep in this 
phase we found evidence for scaling. Furthermore, for given scaling curve, 
increasing N means decreasing K2. Hence, increasing A^ at fixed K2 brings the 
system deeper in the elongated phase. This leads to the conclusion that there 
is nothing wrong with the elongated phase. It decribes very large spacetimes 
which are two dimensional on the scale of r^ but not necessarily at much 
smaller scales. It could be effectively classical at scales much larger than the 
Planck length but much smaller than r^,. 

This reasoning further suggests that A^ and k,2 primarily serve to specify 
the 'shape' of the spacetime. The tuning K2 ~ ^2 i^ apparently not needed for 
obtaining criticality but for obtaining a type of spacetime. The peak in the 
susceptibility of the Regge curvature could be very much a reflection of shape 
dependence. Most importantly, this suggests that the physical properties 
associated with general coordinate invariance will be recovered automatically 
in 4D dynamical triangulation, as in 2D with fixed topologyH. 



^It is good to keep the numbers in perspective: for example in the Standard Model 
/ = 250 GeV and for a Higgs mass m — 100 GeV or less, £ is 15 orders of magnitude 
smaller than the Planck length, or even much smaller. 

■I-A field theory analogue is Zn lattice gauge theory, which for n > 5 has been found to 
posses a Coulomb phase, a whole region in bare parameter space with massless photons; 
see for example ref. pl||. 
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